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Abstract 

We consider a hierarchy of many-particle systems on the line with 
polynomial potentials separable in parabolic coordinates. The first 
non-trivial member of this hierarchy is a generalization of an inte- 
grable case of the Henon-Heiles system. We give a Lax representation 
in terms of 2 x 2 matrices for the whole hierarchy and construct the 
associated linear r-matrix algebra with the r-matrix dependent on 
the dynamical variables. A Yang-Baxter equation of dynamical type 
is proposed. Classical integration in a particular case is carried out 
and quantization of the system is discussed with the help of separa- 
tion variables. This paper was published in the rary issues: Sfb 288 
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1 Introduction 

In the last decade, much work has been carried out in the study of com- 
pletely integrable systems admitting a classical r-matrix Poisson structure 
(c.f. [[[], 0) where the r-matrices depend only on the spectral parameters. 
More recently, interest has developed in the study of completely integrable 
systems where the r-matrices depend also on dynamical variables || £|, [5|, §]. 
It is remarkable that the celebrated Calogero-Moser system, whose com- 
plete integrability was shown a number of years ago (see e.g. f7|), has been 
found only recently to possess a classical r-matrix of this dynamical type 
The development of this ideology came to the series of interesting papers 



In this paper we study another example of a dynamical r-matrix struc- 
ture. It is already well known that it is possible to provide 2x2 Lax operator 
satisfying the standard linear r-matrix algebra for the Hamiltonian system 
of natural form where the potential U is of second degree by the coordinates 
(see e.g. || [15|, [16|, [L7|, |18|). We consider a hierarchy of one-dimensional 
many-particle systems with polynomial potentials of higher degree whose 
complete integrability in the framework of the Lax representation is known, 
but for which the associated r-matrix Poisson structure has not been dis- 
cussed. The systems represent a generalization of the known hierarchy of 
two-particles systems with polynomial potentials separable in parabolic co- 
ordinates, which have the form (see, e.g. 



[N/2] 

v^E^fV)^, (i.i) 



fc=0 



where the positive integer N enumerates the members of the hierarchy. We 
study an integrable many particle generalizations of ( |1 . 1| ) in the framework 
of the Lax representation, written in terms of 2 x 2 matrices, 

tt N \z) = [M N (z),L N (z)}, 

M*> = (wit) %\ «-U°W I ) ■ < L2 » 
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where the functions U(z), V(z), Wn(z), Qn(z) depend rationally in the spec- 
tral parameter z and some constraints are imposed on them. In the first 
nontrivial case (N = 3) the system can be considered as a many-particle 
generalization of a known integrable case of the Henon-Heiles system |H] . 
We emphasize, that the ansatz (|1.2| ) is a natural generalization of the Lax 



representation found in j22] to describe an integrable case of the Henon- 
Heiles system. See also jl7|,|^8| for the link to the su(l, 1)-Gaudin magnet 
which corresponds to a free n-dimensional particle separable in parabolic 
coordinates. 

The polynomial second order spectral problem associated with the Lax 
representation was studied in [23j in the framework of a (K x K) Lax rep- 
resentation, and the recursion relations between different members of a hi- 
erarchy of Lax matrices was also given, but no explicit treatment of any 
associated dynamical system was discussed. The (2 x 2) Lax representation 
we describe appears to us to be more useful to investigate the class of in- 
tegrable systems under consideration and to elucidate the classical r-matrix 
Poisson structure of the system. Such Lax representations were introduced 
in |p4||,|25| and their connection with the restricted coupled KdV flows was 
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(see also [22] for the isomorphism with KdV). But 



studied in 

the r-matrix Poisson structure asosiated with higher degree potentials has 
not been discussed. 

We consider the system within the method of variable separation [29, 
|16| , |17] that permits us to develop the classical theta-functional integration 
theory and to consider the associated quantum problem. The last problem is 
reduced to a set of multiparameter spectral problems which are a confluent 
form of ordinary differential equations of the Fuchsian type. 

The central result of the paper is the description of the Poisson structure 
of the system by a dynamical linear r-matrix algebra defined in V 2 ® V 2 

15,1, 

{ L [ N \x) ® L?\y)} = [d£\x,y),L[ N \x)) - [d™(x,y)zJM(y)], (1.3) 

where ® is the direct product <g> of two matrices, but with the product 
of two matrix elements replaced by their Poisson bracket. The L matrices 
are defined in V 2 <g> V 2 by l[ n \x) = ® I, L ( 2 N \y) = I ® K N \y), 



I is the 2x2 unit matrix, d!yp and d^l 3 are matrices depending both on 
spectral parameters and on dynamical variables through the factor (Qn(x) — 
Qn{v))/(x - y). 
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In the space V 2 ® V 2 ® V 2 , the matrices d\f satisfy the constraints 
(classical Yang-Baxter equation of dynamical type) 

[d[V(x,y),d{f(x,z)} + [d£\x,y),d£\y,z)} + [d£\z,y),d{%\x,z)} + 

+{L{ N \y) ® d[ N 3 \x,z)} + {Lf\z) ® 45V, 2/)} + 

+ [c^\x,y,z),Li N \y)-Li N \z)] =0 (1.4) 

plus cyclic permutations. The matrices gQ^ and c in (|1.4j) involve the dy- 
namical variables through the factor Qn. We follow the now standard but 
somewhat confusing notation that , d\j refer to different matrices de- 
pend on whether the current space is V 2 ®V 2 or V 2 ®V 2 ®V 2 . In V 2 ®V 2 ®V 2 , 
L[ N \x) = L( N \x)®I®id, Lf ] = I®LW(y)®I, L { 3 N) (z) = (z) <g> / ® I . 
The matrix d[^ in V 2 ®V 2 ® V 2 acts like d^ in ith and jth space and as 



/ in the third space. Except for the final term, eqn. ( |1.4| ) is the same as 
that given in ||. The final term in the dynamical Yanng-Baxter equation 
presents as it is emphasised in some new matrix, for which a more general 
ansatz is given in |§. 

A main theme of the paper is to present a new solution of the dynamical 
Yang-Baxter equations (|1.4j) which is associated with a hierarchy of one- 
dimensional dynamical systems. At the same time the r-matrix structure of 
parabolic coordinates is elucidated. 

The main results were announced in the paper |3(| . The classical systems 



separated in ellipsoidal and spheroidal coordinates are concidered in [EHI] 



2 The Hierarchy of Separable Systems 

In this section we describe a hierarchy of completely integrable one-dimensional 
many-particle systems with polynomial potentials. The different members of 
the hierarchy are connected by recurrent relations. We give the Lax rep- 
resentation in terms of (2 x 2) matrices for all the hierarchy, describe the 
associated algebraic curves and present explicit formulae for the integrals of 
motion. 
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2.1 Hierarchy of Hamiltonian Systems with Polyno- 
mial Potentials 

Let us consider the hierarchy of the Hamiltonian systems of n + 1 particles 
defined by the Hamiltonians 

n+l 

H N (p 1 ,...,p n+1 ;q 1 ...,q n+1 ) = -^J), 2 + W w (gi,...,9n+i), (2.1) 

where the potentials Un for each member of the hierarchy are given by the 
recurrence relation 

i n N 

U N = (q n+1 -B)U N _ 1 + -J2Y,(- i y ( lpN- J Ar 2 (2.2) 

4 i=l j=2 

with the first trivial potentials given as 

U = 

U x = -2q n+1 -2B, 

U 2 = -2ql +1 -Wql (2.3) 
z i=i 

At Ai = 0, i — 1, . . . , n, B = the recurrent relations (|2.2|) can be solved as 

[N/2] / n \k 

V N (q u ...,q n+1 ) = - J2 2 1 - 2fc ( JV t) £rf (2-4) 

k=0 \i=l ) 

yielding the principal term of Un(qi, • • • , Qn+i) which is a many-particle gen- 
eralization of ( P-.l| ). The lower degree terms are introduced so as to make 
the many-particle generalization non-degenerate. The expression for the po- 
tentials Un at N > 2 can be written in the form of (N — 2) x (N — 2) 
determinant 





/ 


fN 


9-i 


9o 9i 


• • • 9N-5 


\ 






I'n-i 


-1 


9-i 9o 


■ ■ ■ 9N-6 




U N = (-l)" -1 det 




J'n-2 





-1 9-i 


■ ■ ■ 9N-7 








h 





... 


-1 9-i 








h 







-1 


/ 



(2.5) 
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where 

9-i = q n +i - B, g m = — j — E A Tlh m = 0, . . . 

4 i=i 

2 

fk = Y,U9k-i- 2 ,k = 3,...,N, (2.6) 

1=0 

and the Uq, Ui, IA2 are given by (|2.3|). The first nontrivial potentials are 



U, = -2ql +l - q n+l itli+lit A d + 2Bq 2 n+1 , (2.7) 
1 / ™ \ 2 3 n n / 1 \ 

W 4 = -g (E Qi) - 2^+1 E ^ - 2 <?n+l + E ^ (tfn+l - 2 A ) 



+B ( 4<£ +1 _ 2Bq 2 n+l + q n+1 £ gj ) . (2.8) 

3=1 



The potential ( |2.7| ) is exactly the many-particle generalization of a known 



integrable case of the Henon-Heiles system for which n = 1 (c.f. |32|, |22[). 
Analogously the potential ( |2.8| ) is a many-particle generalization of a "(1 : 
12 : 16)" system known to be separable in parabolic coordinates (see e.g. 



The system with the potential (|2.7| ) possesses the following interesting 
reductions: a) at q n+ i = const, it reduces to the Neuman system, which 
describes the motion of a particle on a sphere in the field of a second order 
potential, and b) at q n+ i = Yh=i $ ^ reduces to an anisotropic oscillator in 



a fourth order potential (see, for instance ||20||) 



q\ - 2 QkQi + Mi = 0, i = 1, • • • , n. (2.9) 



k=l 



The same reduction can be carried out for other members of the hierarchy. 

2.2 The Lax Representation 

We look at a Lax representation of the form ( |1.2j ) with 

71 2 

U(z) = iz-Aq n+1 + AB-J2^-, (2.10) 

i=l Z + Ai 
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V(z) = -\u{z) 



= 2p. +1 + E^- ; (2-11) 

i=l Z + A i 

VMz) = -^(z) + C/(z)Qjv(z), (2.12) 

where Qn(z) is a polynomial of degree N — 2. The ansatz for the function 
t/(z), V(z), W(z) is a generalization of the corresponding ansatz constructed 
by Newell et al. []22] to give the Lax representation for the integrable Henon- 
Heiles system. Here we introduce additional degrees of freedom n > 1, and 
consider higher degrees of the polynomial Q^. See also |17|, |18[ for the link 
of such ansatz to the su(l, 1)-Gaudin magnet which corresponds to a free 
n-dimensional particle separable in parabolic coordinates. 

Proposition 1 The Lax representation ( \L.2j is valid for all the hierarchy of 
Hamiltonian systems ( \2. 1\ ), \2.h\ ) with the polynomial Q n{z) and the function 
Wn{z) given by the formulae 

Q N {z) = zQn^z)-- , (2.13) 

2 oq n+1 

W N (z) = W+(z) + W~(z), (2.14) 

W+(z) = zW+_ 1 -2U N ^ 1 , N = 2,..., (2.15) 

n 2 

W-(z) = £-^r, (2-16) 

i=l z + A i 

where Un-\ is the potential fixing the N — 1-th member of the hierarchy. 

Proof. The Lax representation ( |1.2| ) with the functions ( [2.10| )-( p.l2|) is 
equivalent to the two equations 

W N (z) = V(z)+Q N (z)U(z), (2.17) 
W N (z) = 2Q N (z)V(z). (2.18) 

Substituting ( |2 . 1 0| , p7TTj ) and ( |2.13| ) into (|2.17| ) we obtain after some simplifi- 
cation 

/ \ Pi ( TTl n^, / \ r.dUN-1 Qi dUN-i\ 

Wn[z) = SiTA +z r (z)Q ~- lW " 2 ^r"S ^A % J 



7 



+ *E 



E 



A' 



, 8Un-i ldUN-i^ gf 
+ 2(g n+1 - J B) — + ^ — E 



9g n+ i 2 <9<? n+ i ^ z + Ai 



(2.19) 



Imposing the conditions ( |2.14| ), Q2.15|) , fl2.16| ) we arrive at the following re- 
currence relations for the potentials U^iqi, ■ ■ ■ , q n +i) 



dU 



N 



dqi 
dU N 



1 dU N -i dU N -i . _ 

~q% Ai — - — , i — i, 



2 dq n 



+i 



dqi 



dq 



n+l 



1 " dUn-i m dU N ., 



i=l 



dqi 



dq 



(2.20) 
(2.21) 



n+l 



where we use (2^3) to start the recurrence relations. 

To prove that the equation ( |2.18|) is also satisfied, we proceed by induc- 
tion. First assume that it is valid for N — 1. Then we have 



W N (z) = 2zQ N ^(z)V(z) - (2 Pn+1 + J2 -^L- 

oq n +i \ i= i z + Ai t 

The left hand side of ( [2.22j ) can be rewritten in the form 



(2.22) 



zW+_ 1 (z)-2U N „ 1 +Y, 



Pt 



z + Ai 



zWn-i + 2z}^ — A — 2p n+1 - 



^ z + Ai dqi 
" 8U N -i " Pt &U N 

~ 2 l^Pi^: 2 E 



dq n 



+1 



i=l 



dqi 



r[z + Ai dqi 



Using induction and the equality (|2.20|) completes the proof. 

We can obtain explicit formulae for the functions Wn(z) and Qn{ z ) 



Qn(z) 
W N (z) 



JV-2 



1 dUN^k-i k 



- i: 



Az 



2 fc=0 dq n+ i 

v 1 2 p * 



fc=0 



(2.23) 
(2.24) 
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It easy to see that 

l dW N (x) 

Q N (x) = --j— (2.25) 
For example, for the first nontrivial cases we have 

Q 3 (z) = z + 2q n+1 (2.26) 

n n 2 

W 3 (z) = Az 2 + 4zq n+1 +4Bz + 4q 2 n+1 + J2ql + T,-TT ( 2 - 27 ) 

k=i i=i 2 + A * 

for the many-particle Henon-Heiles system and 

1 n 

Q 4 (*) = z 2 + 2* gri+1 + 3g 2 +1 + -$> 2 -25g n+1 , (2.28) 

8=1 

n 

2 



W 4 (*) = 4z 3 + 4£z 2 + 4z 2 q n+1 + Azq 2 n+l + 2 £ g 

n n 

+ 4g^ +1 + 2 g n+1 E 9i ~ E A d ~ 45 ?n+i 



i=i i=i 



n 2 

+ ErrV ( 2 - 2 9) 

»=1 Z + A « 

for the system with iV = 4. 
2.3 Integrability 

The Lax representation yields the hyperelliptic curve C^ N ' = (w,z), 

Det(L (iV) (z) ~wl) = (2.30) 
generating the integrals of motion H N , F$,i = 1, . . . , n. We have 

w 2 = ]^i(L {N \z)) 2 = V 2 (z) + U(z)W(z) (2.31) 



/From (ggg) and (gig)- (gig) we obtain 
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w 2 = 16^- 2 (^ + 5) 2 + 8#jv + E iV = 3, . . . , (2.32) 

i=l Z + A i 



where 



F$ = 2q* £ ( D ? A'Us , + ^Pn+iPiQi - p- {A* + 4g n+1 - 45) 

k,m=n j2 

+ £ T^4"' z = 1 '---' n (2 - 33) 

k,m=l,k^m m ™ 

with /fcm = 9fcPm — ^mPfc- The coefficients of the curve are the integrals of 
motion, so noting that the coefficient of z"- 1 is 32B and the coefficient of 
z N - 2 is 16B 2 we obtain the expressions for the first potentials ( |2,3| ). Equating 
all the other coefficients of powers of z to zero, i.e. the coefficients of z\ j = 
N — 3, . . . , 1 we reproduce all the hierarchy of integrable potentials in the 
form flOl) . 

Proposition 2 The integrals Hn, F$, i = l,...,n are independent and 
Poisson commute with respect to the standard Poisson bracket, 

{H N , F^} = 0, {Fj>\ F { N k) } = 0, k = 1, . . . , n. (2.34) 

The independence of the integrals is evident. But we postpone the proof 
of ( |2.34|) until §4 where the classical r-matrix structure will be elucidated, 
making the proof trivial. 



3 r— Matrix Representation 

In this section we describe the Poisson structure associated with the Lax 
representation for the hierarchy under consideration. It is found to be a lin- 
ear r-matrix algebra with the r-matrices dependent on dynamical variables 
and satisfying some compatibility conditions which are classical Yang-Baxter 
equations of dynamical type. We also discuss in this section a way to describe 
the system within the standard linear r-matrix algebra by embedding into a 
larger dynamical system. 

Through all the section we use the standard notations for the Pauli ma- 
trices 

a ° =i= {l i)' ai= (? i)' a2= G T)' a3= (J -J- (3 - i} 
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Denote a± = (J\ ±ia 2 - Denote as P the permutation matrix P — \ J2l=o °k ® 
at acting in the product of two spaces V 2 ®V 2 . We denote Pu — P ® I, 
P13 — \ YX=q a k ® I ® °~k and P 23 — I®P the permutation matrices acting in 
the product of three spaces V 2 ® V 2 <8> V 2 . The notation will be used to 
denote an (8x8) matrix acting as a unit matrix in j-th space (i ^ j ^ k). For 
example, let S = cr_£g)cr_, then Si 2 = S®I, S23 = I®S, S13 = cr_Cg>/(g>a-_. In 
addition we introduce the notation ^(x—Xj) = 2Pij/(xi—Xj), s\^\xi,Xj) = 
2a(xi,Xj)Sij, where a(xi,xj) is a scalar function depending on dynamical 
variables. 

3.1 The Classical Poisson Structure 

The following proposition holds: 



Proposition 3 The classical Poisson structure for the hierarchy (\2.1Q, \2.1\ 
2.12) is written in the form 



{L^\x),Lf\y)} = {r{x-y),Lr\x)+Lf\y)]^s{x,y),Lr\x)-Lf\y)l 

(3-2) 

where l[ N \x) = I ® L^ N \x), L% = L,( N \x) (g> I , the matrices r(x — y) and 
S]y(x,y) are given by the formulae 



.(2V) 



r(x — y) 



x-y 



( 1 \ 

10 

10 

V 1 J 



s N (x,y) = 2a N (x,y)S, S = cr_ <8> <x_ 



(3.3) 
(3.4) 



with 



a N (x,y) 



Q N (x)-Q N (y) x N - l -y N - 1 



x-y 

N-3 „,N-k-l 



E 

k=o 



X 



v_ 

x-y 



x-y 
dU 



N-k-l 



dq n 



+1 



(3.5) 
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Proof. Let us rewrite the relation ( |3.2| ) in extended form 



{U(x),U(y)} 
{W N (x),W N (y)} 

{V(x),U(y)} 
{W N (x),U(y)} 
{V(x),W N (y)} 



{V{x),V(y)} = 0, 
4a N (x,y)(V{x)-V{y)), 

-(U(x)-U(y)), 



2 


y 




X 




4 




X 




y 




2 




X 




y 



(3.6) 
(3.7) 

(3.8) 



(V(x)-V(y)), (3.9) 
(W N {x) - W N {y)) - 2a N {x, y)U(x) (3.10) 



with respect to the standard Poisson bracket. The equalities ( |3.6| , p78| , pl9D can 
be proved directly from the definitions Q2.10| , 12.1l| , |2.12|) . Let us prove (|3.7p 
by induction i.e. let ( |3.7| ) be valid for the number N. Then for JV + 1 we 

have 

{W N+1 {x),W N+1 {y)} = {W+ +1 {x), W~(y)} - {W£ +1 (y), W'{x)} 
= {xW+(x) - 2U N , W~(y)} - {yW+(y) - 2U N , W~(x)} 
= x({W+(x), W-(y)} - {W+(y), W~(x)}) 
-y({W+(y), W-(x)} - {W+(x), W~(y)}) 

= x{W+(y), W~(x)} - y{W+(x), W~(y)} - 2{U N , W~{y)} (3.11) 



the inductive hypothesis gives 

{W N+1 (x),W N+1 (y)}-4a N+1 (x,y){V{x) - V(y)) = A 



(3.12) 



with 



2E 



1 



(x + A^iy + Ai) 
dU N 



(2piqi{xQ N {y) - yQ N {x)) 



+ 2(x-y) 



dqi 



- pi x(y + Ai) " Kyj -y(x + Ai)- x 



dqi 



dqi 



(3.13) 



i=l 
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To prove that A = we substitute the expressions fl2.23| ), ( |2.24j ) into A». We 
have 



Aq t xy N -> - yx N -* - i £ (sy* - yx^'"^ 

2 H. oq n +i 



AT-2 



+ 4 J] (^(y + A i )y k -y(x + A i 
dU N 



)x 



k=0 

+ Mx-y) 



dq r , 

dU.N-k-1 
dqi 



dqi 



(3.14) 



^From ( |2.20|) it follows that all Aj = and the derivation of ( |3.7|) is complete. 

Eq. ( ^.1(J[ ) we also prove by induction. Assume that it is valid at N. 
Then 

{V(x),W N+1 (y)} = y{V(x),W N (y)} + (l-y){V(x),W-(y)} 
- 2{V{x),U N } 

2!l ( W £(x)-W+(y)) + ^-(W-(x)-W-(y)) 



x — y " " " x — y 

- 2yU{x)a N {x, y) + (1 - W^j/)} 

- 2{V(x),W w } 

Using the inductive proposition and the formulae ( |2.14| . |2.15] . |2.16| ) we find 

{V(x),W* +1 (y)} 2 
with 



y-a: 



(W^r(x) - W N (y)) + 2a N+1 (x, y)U(x) = A 



A = 2 ^^(^-(x)-^(y)) + (l-y){y(a;),iy-(y)} 
x-y 

- 2W + + 2U{x)Q N {x) -2{V{x)Mn] (3.15) 

The first two terms in (|3.15 ) cancel due to the definitions (p.ll] , pT6[) . To can- 
cel the rest we compute the Poisson bracket directly and use Q2.12 ) . Therefore 
( |3.10| ) is proven. 

The equality ( |3.2|) contains all the information concerning the hierarchy 
of dynamical system under consideration. 
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Let us write the relation ( |3.2| ) in the form 

{4*°(*) ® L?\y)} = [d^\x,y),L[ N \x)] - [d£\x,y)Li N \y)), (3.16) 

with d^j ^ = Tij ~\~ s^jj ^ = S{j — Tij dut i <C j. 
Then@ ' 

{(4*°0«0) fc ® (4%))'} = [d^ l \x,y),Lf\x)} - [d£> k %,y)L?\y)l 

(3.17) 

with 

d£ M (x,y)= (3.18) 

p=0 <j=0 

d^' k,l \x,y) = 

EE(4^(x))^ 1 (4^(y))^ 1 4^(x,y)(4^(x))*(^(»)) ff - 

p=0 <?=0 

As an immediate consequence of ( |3.17| ) and we obtain the proof 

that the conserved quantities H^, F$ are in involution. We have that 

{Tr(L?\x)f,Tr(Li N \y)f} = Tt{(l[ n \x)) 2 ® (4%)) 2 } (3.19) 

Applying to the last equation the equality ( |3.17| ) at k — I — 2 and taking the 
trace we obtain the desired involutivity. 

The Lax representation (|1.2f) also can be recovered from (|3.2| ). We have 

LW(x) = I Jim Tr^^x) <? ^(l/)) 2 }, (3.20) 

where the trace is taken over the second space. Applying the identity ( ft.lTQ 
at k = 1,1 = 2 to ( |3.20| ) we obtain the Lax representation L^ N '(x) = 
[MW(x),LW(x)] with the matrix M^ N \x) given as 

M^ix) = 21im^ 00 Tr 2 4 A °( 2/ )(r(x - y) - s^ N \x, y)). (3.21) 

After the calculation in which we take into the account the asymptotic as- 
sumptions 

lim ^M = 4 hm ^1 = 

y y 

lim^ - U(y)Q N (y)) = 
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we obtain the Lax representation in the form ( |1.2| ). Analogously the Lax 
representation for the higher flows can be obtained. 



3.2 Jacobi Identity 



To prove that the algebra ( |1.3| ) is associative we write the Jacobi identity in 

the V 2 ® V 2 <g> V 2 as 

{L[ N \x) ® {Lf\x) <? 4%)}} + {4*°(*) ? {4^*0 ? 4 V) (?/)}} 
+{Li N \y) ? {4*°(*) <? 4^)}} = (3.22) 

with 4^0*0 = ® / ® J, 4%) = * ® £ (JV) (^) ® 4^0*) = 

I ® J <g) L (Ar) (^). Let us rewrite (|3.22|) in the form 0, 

[4*°0«0, [djJPfoy),^ (*,*)] + [d^\x,y),d^(y,z)} + 
+[d£\z,y),d{¥(x,z)))+ (3.23) 

+ [L[ N \x),{L { 2 N \y) ® - {4*°(*) ? 4? (*,!/)}] + 

+ cyclic permutations = 

We shall show below, that in the cases considered the solutions ( |3.23| ) satisfy 
the three equations 

[d£\x,y),d%\x,z)) + [d%\x,y),d%\y,z)) + [d£\z,y),d£\x,z)) + 

+{Lf\y) ® d%\x,z)} - {Li N \z) ® S»\x,y)} + 

+ [cW (x, y , z) , 4^ (y) - 4^ (*)] = (3.24) 

where c^ N '(x,y, z) is some matrix dependent on dynamical variables. The 



other two equations are obtained from ( |3.24| ) by cyclic permutations. We re- 
mark, that validity of the equations (|3.24 ) with an arbitrary matrix c- N ) (x, y, z) 
is sufficient for the validity of fl3.22| ) and therefore ( |3.24| ) can be interpreted 
as some dynamical classical Yang-Baxter equations. These equations have 
an extra term [c, 4^ — Lj] in comparison with the extended Yang-Baxter 
equations in ||. 

Proposition 4 The following equality is valid for all the member of the hi- 
erarchy of dynamical systems 

{L { 2 N \y)® s 13 (x,z)} - {Li N \z)® s 12 (x,y)} = 2p N (x,y,z)[P 23 ,S 13 + S 12 ] 

d(3 N (x } y } z) (7V) (/v) 



dq, 



n+l 



[s,L^>(y)-L^>(z)} (3.25) 
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with cyclic permutations. In . 23{ ) the matrix s = er_ <S> cr_ <g) a~_ and 
- Q7v(x)(w - 2) +Q A r(w)(2;-x) +Qjv(^)(x-y) 

PN[X,y,Z) — -; — — r 

{x - y)(y - z)[z - x) 
Proof. Let us write the equality (|3.25| ) in the extended form 



{Q(x),Q(y)} 
{V(x),Q(y)} 



{U{x),Q{y)} = 

l da N {x,y) 
4a N {x,y) - - — U(x) 



2 dq 



(3.26) 

(3.27) 
(3.28) 



n+l 



{Q(x),W(y)} + {W(x), Q(y)} = V) (V(x) - V(y)), (3.29) 



dq n 



+1 



{V(z),a N (x,y)} 



4 



-(a N (z,x) - a N (z,y)) 
x-y 

1 U(z) ( dajq{z,y) da^i^z^x) 



2x-y\ dq n . 



dq, 



n+l 



(3.30) 



(157271) is trivial and (gg| ) and (|3~29|) follows immediately from (|3~T0|) , (Q 
and ( |2.25| ). To prove ( |3.30| ) we write it using the properties of H / at(^) and 
V(x) in the form 



E 



Pi 



J x (y + Ai)(x + Ai) 
1 ^ Mi 



(9o;Ar(z, x) daN(z,y) 



dqi 



2^ (2 + + Ai) V dq n+1 

Using the decomposition ( |3.32| ), the equality 



dq 



n+l 



(3.31) 



;r- 



z — £ 



A- h 

z — y 
z-y 



z — X 



z-y 



and the recurrence relation (|2.20|) we find that the equality (p.31|) is valid. 
Note that ( 3.28|) can we rewritten in the form 



{W N (y),a N (z,x}} - {W N (x),a N (z,y)} 
V(x) -V(y) fda N (z,y) da N (z,x) 



(3.32) 



x-y \ dq n+ i dq n+1 J 

All the preceding material can be formulated in the following theorem 
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Theorem 1 Let 

{L[ N \x) ® Li N \y)} = [d[?(x,ylL{ N \x)) - [d™{x,y),lJW{y)) (3.33) 



be the r-matrix algebra constrained by the three conditions on the matrices 
the first of which is written as 

[d£\x,y),d£\x,z)) + [d^\x,y),dS\y,z)] + [d£\z,y),d£\x, z)) + 
+{L[ N \y),d[ N 3 \x,z)} - {Li N \z) ® d£\x,y)} + 

+[ C W(x,y,z) ® 4^(y) - = (3.34) 

and the two other are obtained from ( 3. 34 ) by cyclic permutations. Then 
there exists a solution of the equations ({3. 34\ ) describing the dynamics of the 
hierarchy of completely integrable systems in the n + 1 dimensional phase 
space (pi, gi, . . . , p n+ i, q n +i) given by the formulae 



dif\x,y) 
d)f\x,y) 
c {N) (x,y,z) 
with 



P 



Q N (x) - Q N {y) 



x-y 
P 



S, 



x-y 



(3.35) 



„ + Q^)-Q«(v) Sij ,i< j ,i, j = l,2,3, 



x — y x — y 

d QN{x)(y - z) + Q N {y){z -x) + Q N (z)(x - y) 



dq 



n+l 



Qn{z) 



JV-2 



(x - y)(y - z)(z - x) 
s 3 dU^-k-i k 



k=0 



dq n 



+i 



(3.36) 



and the potentials Un defined as ( \2.b] ). The associated Lax representation 
has the form (l.i) with the matrix elements given 



3.3 Embedding into a Larger System 

In conclusion we point out another direction in the algebraic interpretation 
of the hierarchy of dynamical systems considered. The hierarchy can be 
embedded into a system with more degrees of freedom, which admits the 
standard linear r-matrix algebra (s = 0), 

{L[ N \x) ®Li N \y)} = [r(x-y),L[ N \x)+Li N \y)}. (3.37) 
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To describe this embedding we introduce the following ansatz for the 
functions U(z), V(z), W(z) 



N-2 n 2 

U(z) = E^-ErxT» ( 3 - 38 ) 

4=0 4=1 Z + A * 

N-2 n 

V{z) = E^ + E^TT^ (3-39) 

4=0 i=l Z + 

N-2 n 2 

= 4^- 1 +£wfV + £^-, (3.40) 

4=0 4=1 Z + A i 

where the coefficients Ui, Vl, Wj; depends on the variables q n+ i,p n+ i and in 
addition N — 2 new canonically conjugated variables Qj, Vj, j = 1, . . . N — 
2. The substitution of the ansatz ( |3.38|) - (|3.40|) into ( p.37|) leads to the Lie 
algebra 

{U h Uj} = {V l ,V 1 } = {W l {N \wj N) } = 0,i,j = 0,...,N-2, 
{U h V k } = -2U l+k+1 9(k + l-N-l), 
{V h wi N) } = 2W l+k+1 6(k + l-N-l), 
{w} N \TJ k } = AU l+k+1 9(k + l-N-l),k,l = 0,...,N-2, 
{Vi,wW_ 2 } = 8, i = 0,..., N-2, (3.41) 

where 9(n) = 1 if n > and 9(ri) = if n < 0. The representations of this 
algebra give the algebraic description of the hierarchy of integrable systems 
considered. 

Let us discuss, as an example, the case N = 3 - many-particle Henon- 
Heiles system. We have only the coefficients U ,Ui,V ,Vi, Wq 3 \ in the 
ansatz (|3.38| )- (|3.40|) and the Lie algebra for them is 





= {V ,V 1 } = {W ) 3) M i) } = 0, 


{£>i,Vi} 


= {u 1 ,wl 3) } = {v 1 ,wi 3) } = o, 


{V ,U } 


= -2U 1} 


{wP,u } 


= m, 


{VoM 3) } 


= 2W h 




= {V ,W} 3) } = 8. 



(3.42) 
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The following representation for the algebra ( |3.42| ) can be found 

U = AB - 4g n+1 , Ui = 4, 

V = 2BV l + 2p n+1 -2V l q n+l ,V 1 = 2V l , 



W, 



(3) 



Wi 



(3) 



4g^ +1 + Qi-2P 1 p n+1 + P 1 2 g„+i, 
45 + 4g n+1 -P 1 2 .. 



(3.43) 



The corresponding enlarged dynamical system admits the integrals of 
motion X 3 , TC3, T\,i = 1, . . . , n given by the formulae 



Is 

n 3 



m=l 



1 3 (B - q n+1 ) +V 1 (AJ2 PiQi + Pi E ft' 



i=l 



8=1 



(3.44) 
(3.45) 



*r - 7^3 E a + n 2 (^ Eft 2 - oE ^ft 2 - E ft 2 ) 

8=1 i=l 1=1 8=1 

1 / n n ^ \ 

+2^1 1 ( 5 - ft+i) E^ft _ E ^ft + E ft 2 J ( 3 - 46 ) 

- \ i=l 8=1 i=l / 

where in ( |3.45| ),( p.46| ) H% and i 7 ^ are the integrals of motions of the many- 
particle Henon-Heiles system calculated by the formula ( |2.33|) . We can see, 
that at V\ = 0,X3 = the system is reduced to the many-particle Henon- 
Heiles system. 

The question of the explicit description of the other members of the hi- 
erarchy within this approach will be given elsewhere. 



4 Separation of Variables 

The separation of variables (c.f. | I6| , |TTf ) is understood in the context of the 
given hierarchy of Hamiltonian system as the construction of n + 1 pairs of 
canonical variables 7Tj, //j, % — 1, . . . , n + 1, 

{/ij, /ifc} = {iTi, TT k } = 0, {7Ti, fl k } = 5 ik (4.1) 

and n + 1 functions $j such that 

$, [ N , nj, H N , Fg\ . . . , 4 n) ) = 0, j = 1, 2, . . . , n + 1, (4.2) 
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where Hn, F$ are the integrals of motion in involution. The equations 
are the separation equations. The considered integrable system admits a Lax 
representation in the form of (2 x 2) matrices ( |1.2j ) and we will introduce the 
separation variables 71"$, Hi as 

n = V{m), U(fii) = 0, t = l,...,n + l. (4.3) 

Below we write these formulae explicitly for our system. 

4.1 Parabolic Coordinates 

The set of zeros fij,j = 1, . . . n + 1 of the function U (z) in the Lax represen- 



tation (|1.2f) defines the parabolic coordinates given by the formulae P9| , \L7 



n n+1 

Qn+1 = ^2 A + B + Hi, 
i=l i=l 

2 . TYjtKVj + An) 1 ., ^ 1 (a a\ 

q m = -4 J — , m = l,...,n, if n>l (4.4) 

and 

qf = -A(fx 1 + A)(fi 2 + A),A 1 = A, if n=l. (4.5) 
Let us denote by 7r m , 

vr m = V(fii) = (i m II ~ ~~T > m = 1 > • • • > n + L ( 4 - 6 ) 

i=l,. ..,»+! 



Proposition 5 TTie coordinates Hh^i given by fal-$) > fii^ ) are canonically 
conjugated. 

Proof Let us list the commutation relations betveen U(z) and V(z), 



{U(x),U(v)} = {V(x),V(y)} = 0, (4.7) 

(4.8) 
(4.9) 



{V{x),U{y)} = -?—(U(x)-U(y)), (4.8) 
y-x 
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The equalities {fii, fij} = 0, i,j = 1, . . . ,n + 1 follows from ( (4.7[) . To derive 



the equality {/Xj, 7Tj} = Sij,i,j = 1, . . . , n + 1 we substitute x = /ij in (|4.8|) , 
obtaining thus 

{ irjt U(y)} = --^—U(y)), 

y 

which together with the equation 

= {tTj, [/(//;)} = {nj,U(y)} \ y= ^ +U'(fi i ){n j , //J, 

gives 

fa, //J = -__fa,[%)} | y=Mi = 
Equalities {77,, ttj} = 0,i,j = 1, . . . ,n+ 1 can be verified by the similar way: 

-fa, 71^} = 

{v(^-).^(y)}l«= w +V(ik){t*i,vM} 

V\ N ){V( N )^ 3 } +V'(^){^V(^)} = 0. 

The separation of variable equations have the form 

?r? = w 2 (fMi), i = l,...,n + l, (4.10) 

where the function w 2 (z) is given by ( |2.32| ). Our use below of the separation 
equations is two-fold - to integrate the equations of motion in terms of theta 
functions and to quantize the system. 

4.2 Theta Functional Integration of the Many Particle 
Henon— Heiles System 



The curve ( |2.32| ) has genus g = n+ . So only in the cases N = 3, 4 does 



the number of degrees of freedom coincide with the genus. We consider here 
only the case N = 3 (the Many- Particle Henon- Heiles System) for which the 
curve (|2.32| ) has a branching point at infinity. With this in mind we reduce 



( |4 . 1 U ) to the Jacobi inversion problem 



"+1 mu ,,k, 



2_, / ~~T\ =C k , k = 0,...,n 



i=i 

n+1 



E/"*=«K ("I) 
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where C n+ \ = (y, z) is a nonsingular hyperelliptic curve of genus n + 1 given 
by the formula 



V 



w 



n+l 

2 11(2 + A 

m=l 



2n+3 

fe=i 



*(<&)), ft + Qi ( 4 - 12 ) 

, ■ ■ Q2n+3,Q2n+4, z{Q2u+a) = OO 



with the function w 2 given in ( |2.32|) and Qi 
being the branching points. 

Let (a, b) = (ai, . . . , a n+ \] b\, . . . , b n+ x) be a canonical base of the first 
homology group Hi(C n+ \, Z) with the intersection matrix, 



I 



n+l 




and let v = (v x , . . -,v n+1 ), 

n+l 

V 3 = CjkWg-k+l, 
fc=l 



y 



,n + l 



(4.13) 



be the basis of holomorphic differentials normalized in such a way, that the 
period matrix of v n+ i, . . . , t> n+ i with respect to (a, b) has the form I1 T = 
(l n+ i;r), where l n+ i is the (n + 1) x (n + 1) unit matrix, and the (n + 
1) x (n + l) matrix r belongs to the Siegel upper half-space of the degree 
n + l, «S n+ i = {t I r = t t ,Ihit > 0}. Let D = (A,B) be the divisor, 
where A and B are positive divisors of degree g; let us denote the Jacobian 
variety of the curve C n+ i by J(C n+ \) = C n+1 /(l n+ i, r) and the mapping 
(Abel map) which establishes the correspondence between the point e = 
(Xa v ) g J(Cn+i) and divisors V G C n+i by V => J(C n+1 ). We write 
e = \(e\e") , (e', e") e C n+1 , where [e] = is the characteristic 
of the point. If e'^s" is or 1, the characteristics [e] are the characteristics 
of half-periods. 

We identify the branching points Qj, j = 1, . . . , 2n + 4 with the char- 
acteristic of half-periods e J(C n+ \) and fix for definiteness the basis of 
homologies (a, b) as follows. 



e-i 

e 2n+3 



0...0 
0...0 



e 2 



00.. .0 
10...0 



e 3 



1 10. ..0 
100.. .0 



00...010...0 
1 1...100...0 



10...0 

1 



^2fc+2 
e 2n+4 = 



00...010...0 
1 1...110...0 



10...0 
00.. .0 



(4.14) 
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It is naturally to chose the a-cycles to be homological to the trajectories of fi 
variables under the evolution of the system . The vector of Riemann constant 
K in the basis QOg ) is K = [^Xiki 

The Riemann theta function 0[e](z|r) with the characteristic [e] = 
is determined on C n+1 x S n+ i as 



Q[s}( 



z r 



5' 5' 

= ex P 7rv /Z l{((m+ -)r, (m + -)) 

mez»+i z z 

+ 2((m+ £ -),z+ £ -)}, 



(4.15) 



where (•, •) means the Euclidean scalar product. Let us consider the curve 
with the ordering of the branching points 



Q 2 k = A k , k = l,...,n. 
Using known hyperelliptic theta formulae 



3=1 

where k — 1, 



?n+l(*) 



©(C 4 v + (/S + . . . + J& +1 )v + K; r) 
n and ft^ are the constants we obtain 
2/nA 02 [ e 2"+4 + K ]Oo; r)e 2 [e 2fc+2 + K](ut + uj ;t) 



(4.16) 



(4.17) 



«(0) 



© 2 [e 2 fc+2 + K] (w ; r) 9 2 [e 2n+4 + K] (wi + uj ] r) 

n+l 



7=1 i=l Jai 



d 2 



— ln9[e 2n+4 + K](atf + a; ;r) 



(4.18) 



where the "winding vector" uj is expressed in terms of b-periods of the differ- 
ential of the second kind Q with zero a-periods and a unique pole at infinity 
which is expanded as Q(z) = d(/( 2 + o(l)d(, z = l/( 2 = oo as 



16ni 



n,...,f n 



and ioq G C 9+1 is a constant. 
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The set of angles 9{ = Uit + u>iQ, i = 1, . . . , n + 1 and the n + 1 a-periods 
of Abelian differentials 

Jk— f KkdfJ-k = f wdz, k — 1, . . . , n + 1, (4.19) 

constitute the "action-angle" variables J^, 6^ for the system. The standard 
equality 

n+l 

# 3 = £^ (4-20) 

follows from the Riemann bilinear identity. 

The theta functional formulae for the evolution of the system can also be 
given for the case N = 4, where the curve ( j2.32| ) has no branching points at 



infinity, but its genus coincides with the number of degrees of freedom (see, 
e.g. [34]). The algebro-geometric integration of the systems corresponding 



to the higher members of hierarchy (N > 4) requires special consideration. 
4.3 Quantization 

The separation of variables has a direct quantum counterpart |$5|, @. To 
pass to quantum mechanics we change the variables 7Tj,//j to operators and 
the Poisson brackets (14. II) to the commutators 



[/ij,/i fe ] = [n^TTk] = 0, [vr^/ifc] = -i5 ik (4.21) 

Suppose that the common spectrum of /ij is simple and the momenta 7Tj are 
realized as the differentials tcj = —i-^-. The separation equations ( |4.10| ) 
become the operator equations, where the noncommuting operators are as- 
sumed to be ordered precisely in the order as those listed in (|4.2|) , that is 
7Tj, yUj, Hjsf, Fff, . . . , F^. Let • • • , fin+x) De a common eigenfunction of 

the quantum integrals of motion: 

H N V = A n+1 ^, F$if! = A f tf , i = 1, . . . , n. (4.22) 

Then the operator separation equations lead to the set of differential equa- 
tions 

d 

^(-t—, fit, H N , f£ ] , . . . , F^ivi, . . . , Hn+i) = 0, j = l,...,n + l, 



dfjLi 

(4.23) 
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which allows the separation of variables 

n+l 

#0i,..., /wi) = IlV'i(^)- (4.24) 

i=i 

The original multidimensional spectral problem is therefore reduced to the 
set of one-dimensional multiparametric spectral problems which have the 
following form in the context of the problem under consideration 

/ d 2 n A \ 

+ 16x N - 2 (x + Bf + 8A n+1 + —^r A i> • • • > A «+i) = 



(4.25) 

with the spectral parameters Ai, . . . , A n+ i The problem ( [4.25] ) must be solved 
on the n+l different intervals - "permitted zones": x G [Qi,Q 2 ], ■ ■ ■ , 
[<?2n+i, Q2n+2] for the variable x. 
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